Roton-maxon spectrum and stability of trapped dipolar condensates 
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We find that pancake dipolar condensates can exhibit a roton-maxon character of the excitation 
spectrum, so far only observed in superfluid helium. We also obtain a condition for the dynamical 
stability of these condensates. The spectrum and the border of instability are tunable by varying 
the particle density and/or the confining potential. This opens wide possibilities for manipulating 
the superfluid properties of dipolar condensates. 



Recent progress in cooling and trapping of polar 
molecules [l|, |2| opens fascinating prospects for achiev- 
ing quantum degeneracy in trapped gases of dipolar par- 
ticles. Being electrically or magnetically polarized, po- 
lar molecules interact with each other via long-range 
anisotropic dipole-dipole forces. This makes the prop- 
erties of such dipolar gases drastically different from the 
properties of commonly studied atomic cold gases, where 
the interparticle interaction is short-range. Other candi- 
dates to form a dipolar gas are atoms with large mag- 
netic moments |i3-Lll- and atoms with dc-field- or light- 
induced electric dipole moments 0, Q> @ • 

The dipole-dipole interaction is responsible for a va- 
riety of novel phenomena in ultracold dipolar systems. 
The energy independence of the dipole-dipole scattering 
amplitude for any orbital angular momenta provides re- 
alistic possibilities for achieving a superfluid BCS tran- 
sition in single-component dipolar Fermi gases (see 
and refs. therein). Dipolar bosons in optical lattices 
have been shown to provide a highly controllable envi- 
ronment for engineering various quantum phases [lfj . In 
addition to superfluid and Mott-Insulator ones, recently 
observed in Munich experiments with bosonic atoms, 
the long-range dipole-dipole potential provides supersolid 
and checker-board phases. Dipole-dipole interactions are 
also responsible for spontaneous polarization and spin 
waves in spinor condensates in optical lattices |l2| . and 
may lead to self-bound structures in the field of a travel- 
ing wave jl-'i| . Recently, dipolar particles have been con- 
sidered as promising candidates for the implementation 
of fast and robust quantum-computing schemes 0, ITij . 

The long-range and anisotropic (partially attractive) 
character of dipole-dipole forces ensures a strong depen- 
dence of the stability of trapped dipolar Bose-Einstein 
condensates on the trapping geometry Q- For cylindri- 
cal traps with the aspect ratio l z /l p > 1* = 0.43, a purely 
dipolar condensate is dynamically unstable if the number 
of particles N exceeds a critical value. A detailed study 
of the excitation modes for this geometry is contained in 
Ref. ,15|- It has also been argued in Ref. Q that in 
pancake traps with l z /l p < Z* the ground state solution 
is expected for any N. 



In this Letter we analyze the nature of excitations and 
instability of pancake-shaped dipolar condensates. For 
this purpose, we consider the physically transparent case 
of an infinite pancake trap, with the dipoles perpendic- 
ular to the trap plane. For the maximum condensate 
density uq — > oo, the dynamical stability requires the 
presence of a sufficiently large short-range repulsion, in 
addition to the dipole-dipole interaction. Otherwise, if 
no exceeds a critical value n c , excitations with certain 
(large) in-plane momenta q become unstable. At densi- 
ties n < n c the excitation spectrum has a roton-maxon 
form (see Fig.l) similar to that in superfluid helium. 

The roton-maxon spectrum for helium has been sug- 
gested by Landau 0], and later Feynman related 
the excitation energy to the static structure factor of the 
liquid. The roton minimum originates from the fact that 
at intermediate momenta one has a local structure pro- 
duced by the tendency of atoms to stay apart. The pan- 
cake dipolar condensate is the first example of a weakly 
interacting gas, where the spectrum has a roton-maxon 
form. In an infinite pancake trap this spectrum allows a 
transparent physical interpretation. 

We emphasize that the roton-maxon spectrum finds 
its origin in the momentum dependence of the interpar- 
ticle interaction. For in-plane momenta q much smaller 
than the inverse size L of the condensate in the confined 
direction, excitations have two-dimensional (2D) char- 
acter. Hence, as the dipoles are perpendicular to the 
plane of the trap, particles efficiently repel each other 
and the in-plane excitations are phonons. For q 1/L, 
excitations acquire 3D character and the interparticle re- 
pulsion is reduced. This decreases the excitation energy 
under an increase of q. The energy reaches a minimum 
and then starts to grow as the excitations continuously 
enter the single-particle regime. The minimum energy 
is zero for the maximum density no equal to a critical 
value. At higher densities, excitations with momenta 
q in the vicinity of this minimum become unstable and 
the condensate collapses. Below we find the condensate 
wave function, excitation spectrum, and the conditions 
for both rotonization and instability. 

We consider a condensate of dipolar particles harmon- 
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ically confined in the direction of the dipoles (z) and uni- 
form in two other directions (p = {x 7 y}). The dynamics 
of the condensate wave function ip(r, t) in this infinite 
pancake trap is described by the time-dependent Gross- 
Pitaevskii (GP) equation (see and refs. therein): 



f A+> 2 +#(r,t)| 2 
2m 2 



dr'V d (r-r')\^(r',t)\ 2 }^(r,t), 



(1) 



where w is the confinement frequency, m is the parti- 
cle mass, and d the dipole moment. The wave function 
ip(r, t) is normalized to the total number of particles. The 
third term in the rhs of Eq.© corresponds to the mean 
field of short-range (Van der Waals) forces, and the last 
term to the mean field of the dipole-dipole interaction. 
The coupling constant for the short-range interaction is 
g, and Vd(r) = (1 — 3cos 2 #)/r 3 is the potential of the 
dipole-dipole interaction, with 9 being the angle between 
the vector r and the direction of the dipoles (z). 

The ground state wave function is independent of 
the in-plane coordinate p and can be written as 
tpo(z) exp (— ijit), where \x is the chemical potential. 
Then, integrating over dp' in the dipole-dipole term of 
Eq.©, we obtain a one-dimensional equation similar to 
the GP equation for short-range interactions: 

~ :L A + Y wV + (3 + 9d ^(z)-(^ M*) = °> 

(2) 

where g d = 87rd 2 /3. We will discuss the case of (g+gd) > 
0, where the chemical potential \x is always positive. For 
/i ^> Tiuu the condensate presents a Thomas-Fermi (TF) 
density profile in the confined direction: V'oC 2 ) = n o(l — 
z 2 / L 2 ), with riQ = n/(g+gd) being the maximum density, 
and L = (2u/mcu 2 ) 1 / 2 the TF size. 

Linearizing Eq.© around the ground state solution 
ipo(z) we obtain the Bogoliubov-de Gennes (BdG) equa- 
tions for the excitations. Those are characterized by the 
momentum q of the in-plane free motion and by an in- 
teger quantum number (J > 0) related to the motion in 
the z direction. The excitation wave functions take the 
form f±(z) exp(iqp), where f± = u ± v, and u, v are the 
Bogoliubov {u, v} functions. Then the BdG equations 
read: 



ef- 
e/4 



2m 



' dz 2 



Hkinf- + H int [f-], 



Hkinf+, (3) 

(4) 



where Uhin is the sum of kinetic energy operators, and 

R int [f-\ = 2(g d + g)f_(z)^ 2 (z) - (3/2)g d qM^) X 
x / dz' f^(z')Mz')^v(-q\z-z'\). (5) 



For each j we get the excitation energy as a function of 
q. We will be mostly interested in the lowest excitation 
branch eo (q) for which the confined motion is not excited 
in the limit q — > 0. 

The second term in the rhs of Eq.lJSJ originates from 
the non-local character of the dipole-dipole interaction 
and gives rise to the momentum dependence of an effec- 
tive coupling strength. In the limit of low in-plane mo- 
menta qL <C 1, this term can be omitted. In this case, 
excitations of the lowest branch are essentially 2D and 
the effective coupling strength corresponds to repulsion. 
Equations and (J2J become identical to the BdG equa- 
tions for the excitations of a trapped condensate with a 
short-range interaction characterized by a coupling con- 
stant (g + gd) > 0. In the TF regime for the confined 
motion, the spectrum of low- energy excitations for this 
case has been found by Stringari [18j . The lowest branch 
represents phonons propagating in the x, y-plane. The 
dispersion law and the sound velocity c s are given by 



e (q) = hc s q; 



(2/x/3m) 1/2 . 



(6) 



For qL ^> 1, the excitations become 3D and the effec- 
tive coupling strength decreases. The interaction term is 
then reduced to Hi nt \f-\ = (2g — gd)ipo(z)f-(z) as in the 
integrand of Eq. J^J one can put z' — z in the arguments 
of /_ and ifjQ. In this case, Eqs. J2J and (@J are similar 
to the BdG equations for the excitations of a condensate 
with short-range interactions characterized by a coupling 
constant (2g — gd). If the parameter j3 — gj gd > 1/2, this 
coupling constant is positive and one has excitation en- 
ergies which are real and positive for any momentum q 
and condensate density no. For (3 < 1/2, the coupling 
constant is negative and one easily shows that at suffi- 
ciently large density the condensate becomes unstable. 
For collective excitations in the TF regime at no — > oo, 
kinetic energy terms in Eq. Q can be omitted and it reads 
ef + = (2g — gd)no(z)f-. Then, rescaling the excitation 
energies as e 2 = e 2 (2g — g d )/(g + g d ), Eqs. © and (3J give 
the eigenmode equation e 2 /+ = 2(g + gd)no(z)Hki n f+ for 
positive excitation energies e in the case of short-range 
repulsive interaction with the coupling constant (g + gd)- 
Thus, for j3 < 1/2 we obtain e 2 < and imaginary e, 
which indicates dynamical instability of the condensate 
with regard to these high-momentum excitations. 

We thus see that the most interesting behavior of the 
excitation spectrum in the TF regime is expected for 
qL ^> 1 and (3 close to the critical value 1/2. In our 
analytical analysis we first reduce Eqs. © and J2J to 
the equation for the function W defined by the relation 
= Wyl — x 2 , where x = z/L. Expressing the function 
/_ through W from Eq.©, we substitute it into Eq.© 
and integrate straightforwardly over dz' in Hi nt [f-\ as 
the main contribution to the integral comes from a nar- 
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row range of distances 



l/q. This yields 



x 2 ) 



d 2 W 
dx 2 
2/3-1 

i + p 



- 1 



2(1 + 0) 



HE 9 (1 - x 2 ) - 



dW 

dx 

3fiV 



2(1 + /?) 



*2 2 i 

n id + 



W = 0, (7) 



where E q — h q /2m. Here we omitted terms of the 



order of E q % lo < / n and h w / n , since they are small 
compared to either Ti 2 Lo 2 or E 2 . 

For each mode of the confined motion (each quantum 
number j), the solution of Eq.Q can be written as series 
of expansion in Gegenbauer polynomials C^(x), where 
A = (4 + /3)/2(l + /3), and n > is an integer. The cou- 
pling between polynomials of different power is provided 
by the term proportional to (2/3 — 1)(1 — x 2 )W. For the 
critical value (3=1/2 the coupling is absent, and we then 
obtain Wj oc Cj{x). The dispersion law is characterized 
by a plateau (see Fig. la), and for the j-th branch of the 
spectrum it is given by 



^(?)=£ ? 2 +fiV(l+i(j'+3)/2); (3 



1/2, gi» 1. (8) 



For (3 ^ 1/2, assuming that the coupling term 
fiE q \2/3 — 1|/(1 +0) % h 2 uj 2 and it does not significantly 
modify the eigenfunctions, we can still confine ourselves 
to the perturbative approach. Then, as the polynomials 
Cj are orthogonal with the weight (1 — x 2 ) A_1 / 2 , for the 
lowest branch of the spectrum we obtain 



e 2 (l) = El 



(2/3-l)(5 + 2/3) 
3(l + /3)(2 + /3) 



E q ^+h 2 uj 2 ; qL > 1. (9) 



From Eq. @ one sees two types of behavior of the spec- 
trum. For j3 > 1/2 the excitation energy monotonously 
increases with q (see Fig. lb). If j3 < 1/2, then the dis- 
persion law is characterized by the presence of a min- 
imum. Since in the limit of qL <C 1 the energy eo grows 
with q, the existence of this minimum indicates that the 
spectrum as a whole should have a roton-maxon charac- 
ter (see Fig. lb). This behavior is known from the physics 
of liquid helium. As discussed above, in our case it is re- 
lated to the reduction in the coupling strength with an 
increase of momentum, resulting from the transformation 
of the character of excitations from 2D to 3D. 

As follows from Eq.@ for (3 close to 1/2, the roton 
minimum is located at q = (lG/x^/lSTia;) 1 / 2 !/^, where 
S = 1/2 — P, and lo — (h/muj) 1 ^ 2 is the harmonic oscilla- 
tor length for the confined motion. The excitation energy 
at this point is e m i n = [h 2 uj 2 — (S/^/15) 2 ] 1 / 2 . An increase 
of the density (chemical potential) or S makes the roton 
minimum deeper. For fiS/hcu = 15/8 the minimum en- 
ergy reaches zero at q — V2/Iq. At larger values of /iS/hca 
one gets imaginary excitation energies for q ~ l/^o, and 
the condensate becomes unstable. 

We have then found the excitation spectrum numeri- 
cally from Eqs. and (0J for various values of (3 and 



H/Tiw. The results for the TF regime and (3 close to 1/2 
are presented in Fig.l, where one sees a good agreement 
between the numerics and analytics. The discrepancy 
is mainly due to the neglect of the border effects and 
some of the kinetic energy terms when obtaining Eq. Q 
from Eqs. J3J) and iQJ). A similar behavior of the spec- 
trum is observed for non-TF condensates. In this case, 
due to a large kinetic energy in the confined direction, 
the stability of the condensate does not require as strong 
a short-range repulsive coupling strength as in the TF 
regime. Accordingly, the rotonization of the spectrum 
and the instability appear at smaller values of (3. These 
critical (3 have been calculated numerically as functions 
of h/Tiuj and are shown in Fig. 2. 

The dipolar condensate is the first example of a weakly 
interacting gas offering a possibility of obtaining a roton- 
maxon dispersion, up to now only observed in the rela- 
tively more complicated physics of liquid He. In contrast 
to the helium case, the rotonization in dipolar conden- 
sates is tunable. By varying the density, the frequency of 
the tight confinement, and the short-range coupling one 
can manipulate and control the spectrum, making the 
roton minimum deeper or shallower. One can also elimi- 
nate it completely and get the Bogoliubov-type spectrum 
or, on the opposite, reach the point of instability. 

The instability of dipolar condensates with regard to 
short-wave excitations, is fundamentally different from 
the well-known instability of condensates with attractive 
short-range interaction (negative scattering length). In 
the latter case the chemical potential is negative and the 
ground state does not exist. The unstable excitations 
are long-wave and an infinitely large cloud undergoes lo- 
cal collapses. For the dipolar BEC the chemical potential 
is positive and the instability is related to the momen- 
tum dependence of an effective coupling strength. The 
unstable excitations become the ones with high momenta 
at which the coupling is attractive. The existence of the 
roton minimum at a given (3 < 1/2 for fx/huj just below 
the point of instability, is likely to indicate that there is 
a new ground state in the region of the condensate insta- 
bility. The presence and character of this state will be a 
subject of our future studies. 

The presence of the roton minimum in the excitation 
spectrum can be observed in Bragg-spectroscopy experi- 
ments as those of Steinhauer et al. [ljj, or in the MIT- 
type of measurement of the critical velocity for super- 
fluidity [2(|. According to the Landau criterion |2l| . 
the critical velocity v c is equal to the minimum value of 
eo(q)/q, and the presence of the roton minimum strongly 
reduces v c . Even in the absence of rotonization, a de- 
crease in the slope of the dispersion curve at large mo- 
menta leads to a significant reduction of the critical ve- 
locity. 

In conclusion, we have found that pancake dipolar con- 
densates can exhibit a roton-maxon character of the ex- 
citation spectrum. The presence, position, and depth 
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FIG. 1: Dispersion law eo(g) for various values of j3 and [i/hui: 
(a) P = 1/2, fj,/hu> = 343; (b) /3 = 0.53, fi/hiv = 46 (upper 
curve) and /3 = 0.47, [i/fiu) = 54 (lower curve). The solid 
curves show the numerical results, and the dotted curves the 
result of Eq. ©. 
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FIG. 2: Critical values of (3 for the rotonization (filled squares) 
and for the instability (hollow squares) versus fi/hui. 

of the roton minimum are tunable by varying the den- 
sity, confining potential, and the short-range coupling 
strength. This opens new handles on manipulations of 
superfluid properties of trapped condensates. 
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Note added: After this work has been completed, we 
become aware of the recent preprint by O'Dell et al. 
where a roton-maxon spectrum has been obtained numer- 
ically for some particular cases in elongated condensates 
with laser-induced dipole-dipole interatomic interactions. 
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